The correction-to-scaling exponent in dilute systems 
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The leading correction-to-scaling exponent lo for the three-dimensional dilute Ising model is calcu- 
lated in the framework of the field theoretic renormalization group approach. Both in the minimal 
subtraction scheme as well as in the massive field theory (resummed four loop expansion) excellent 
agreement with recent Monte Carlo calculations [Ballesteros H G, et al Phys. Rev. B 58, 2740 
(1998)] is achieved. The expression of lo as series in a -^/e-expansion up to 0(e 2 ) does not allow a 
reliable estimate for d = 3. 
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From renormalization group (RG) theory one knows that 
in the asymptotic region the values of the critical expo- 
nents are universal and scaling laws between them hold. 
There the couplings of the model Hamiltonian describing 
the critical system have reached their fixed point val- 
ues. In the nonasymptotic region deviations from the 
fixed point values are present. They die out according 
to a universal power law governed by the correction-to- 
scaling exponent lo. E.g. for the zero field susceptibility 
the approach from above to the critical temperature T c 
is characterized by the so-called Wegner expansion!!! 

x c~ r r-T (i + y x t"I v + r 2 T 2 "/" + ...), (l) 

where r = (T — T c ) /T c and the Ti are the non- universal 
amplitudes. 7 and v are the asymptotic values of the sus- 
ceptibility and correlation length critical exponents. The 
smaller the exponent to, the larger is the region where cor- 
rections to the asymptotic power laws have to be taken 
into account. Being even further away from the fixed 
point it is necessary to consider the complete non linear 
crossover functions. 

The implication of quenched dilution on the critical be- 
havior is a long-standing problem attracting theoreti- 
cal, experimental and numerical efforts. In the 3e?-Ising 
model quenched disorder changes the asymptotic criti- 
cal exponents compared to the pure onesm In principle 
this statement should hold for arbitrary weak dilution. 
But in order to observe this change one should approach 
the critical point close enough. The width of this region 
turns out to be dilution dependent. 

In particular Monte Carlo (MC) calculations of the crit- 
ical exponents in the dilute 3d-Ising model are more dif- 
ficult to perform than for the pare model since they 
need much larger sizes of latticesa Even then the ex- 



ponents were found to be non-universal and varyingjpon- 
tinously with dilution, i.e. they were effective onestl. It 
became clear that a correction-to-scaling analysis is un- 
avoidable and indeed universal exponents were foundO. 
Without it one still obtains concentration dependent ef- 
fective exponents!! 

The value of the correction-to-scaling exponent u> found 
in MC calculations from an analysis invoking the first 
correction term in (|l|) turned out to beO 

lu = 0.37 ±0.06. (2) 

Thus it is almost half as large as its corresponding value 
in the pure model (see Table |) and this smallness of lo 
in the dilute case explains its importance for an analysis 
of the asymptotic critical behavior. It is therefore highly 
desirable to have an independent quantitative theoreti- 
cal prediction for the value of the correction-to-scaling 
exponent in the dilute system. 

In theoretical calculations the value of lo found by scal- 
ing field RGB is lo = 0.42. So far field theoretical RG 
studies mainly concentrated on the asymptotic values of 
the leading exponents. Correction-to-scaling exponents 



Method 


dilute 


pure 




scaling field 
e expansion 
massive RG, d = 3 
min. sub. RG, d = 3 

MC 


0.42K 

see text 
0.372 ± 0.005 
0.390 ± 0.04 

0.37±0.0efl 


0.870 

0.814 ±0.018 
0.799 ±0.011 
0.791 ± M36 
0.8 ± O.lEAj 
0.8-0.85SL 
0.87±0.09E3 
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TABLE I. Values of correction-to-scaling exponent lo as 
obtained from different methods in dilute and pure 3d-Ising 
models. 
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have been calculated within massive RG in two loop ap- 
proximation in Ref.EI (u) — 0.450) and within the minimal 
subtraction scheme in three loop approximation in Ref.t2l 
(lu = 0.366). Here, we improve this value in the massive 
RG scheme up to four loop order with the result 



u = 0.372 ±0.005 



(3) 



in excellent agreement with (g) . In the minimal subtrac- 
tion scheme we obtain lu = 0.390 ± 0.04 remaining with 
in the bandwidth of MC accuracy. 

The critical behavior of the quenched weakly dilute Ising 
model in the Euclidian space of d = 4 — e dimensions is 
governed by a Hamiltonian with two couplingsEj: 



ml4> 2 a _ 



n 

n 



(4) 



a=l 



in replica limit n — ► 0. Here <p a are the components of 
order parameter; u Q > 0,v Q > are bare couplings; 
is bare mass. 
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We describe the long-distance properties of the model 
(Q) in the vicinity of the phase transition point using a 
field-theoretical RG approach. The results presented in 
this paper are obtained on the basis of two different RG 
schemes: the normalization conditions of massive renor- 
malized theory .at fixedl!3 d = 3 and the minimal sub- 
traction schemeE^L-|The last approach allows both fixed 
d = 3 calculationsli-3 as well as an e-expansion. 

In the RG method the change of the couplings u and v 
under renormalization is described by two /3-functions 



dv 

<W ' 



(•5) 



where /i corresponds to the mass in the massive field the- 
ory approach and to the scale parameter in the minimal 
subtraction scheme. The subscript in (Q) indicates that 
the derivatives are taken at constant unrenormalized pa- 
rameters. The /3-functions differ for different RG schemes 
and in consequence the fixed point coordinates u* , v* , de- 
fined by the simultanious zeros of both /3-functions, are 
scheme dependent. The asymptotic critical exponents as 
well as the correction-to-scaling exponent do not depend 
on the RG scheme and take universal values. 

The correction-to-scaling exponent lu is defined by the 
smallest eigenvalue of the matrix of derivatives of the 
/3-functions 



80 u d@u 
du dv 

du dv 



(6) 



taken at the stable fixed point. For the stable fixed point 
both eigenvalues of this matrix have a positive real part. 

Our results for the correction-to-scaling exponent are 
based on the known high order expansions for the func- 
tions P u and (3 V . In the massive scheme they are known 
in four loop approximatiortJ. In the minimal subtrac- 
tion scheme one can obtain these functions in five loop 
approximation in the replica limit from those of a cu- 
bic modell3. In the limiting case of the pure model only 
the coupling u is present. The corresponding /3-function 
results from putting v = in (3 u (u, v) and the correction- 
to-scaling exponent is simply the derivative d/3 u (u, 0) /du 
taken at the stable fixed point u* . Note that for the pure 
model the /3-functions in the massive scheme are known 
in six loop approximatiorO and the five loop results fpc 
the RG functions in the minimal subtraction schemet^l 
agree with those recovered from Refill 

It is known that the series obtained in the perturbational 
RG approach arerajt best asymptotic (for the dilute model 
see however Refill). An appropriate resummation pro- 
cedure has to be applied to the /3 functions in order to 
obtain reliable information. The choice of the resum- 
mation procedure depends on the information about the 
high order behavior of the expansion series. This infor- 
mation is not available for the case of the /3-functions (||) . 
In this situation the most appropriate way to proceed, is 
to use the Pade-Barel resummationE-3 generalized for the 
two variable caseE3. 

The steps which we follow in the calculation of the 
correction-to-scaling exponent lu are the following: First 
the /3-functions (||) are resummed and the system of equa- 
tions for the fixed points, f3 u {u* ,v*) = 0, @ v 0u*,v*) = 0, 
is solved. Then the matrix of derivatives (o) is calcu- 
lated for the resummed /3-functions. The stability of the 
fixed points is checked. The fixed point with both u* ^ 
and v* 7^ is the stable one at d = 3 and the smallest 
eigenvalue gives the desired correction-to-scaling expo- 
nent. Note that the eigenvalues might be complex, in 
this case both have the same positive real part defining 
to. 

In Fig. we present our results for the exponent lu 
obtained in successive orders of perturbation theory in 
number of loops. To perform the resummation the Borel 
transforms of the truncated Zth order perturbation theory 
expansion for the /3-functions were presented in the fopa 
of [(I — 1)/1] rational approximants of two variablesE3. 
This form of rational approximants appeared to give the 
most reliable results. The four loop results for the ex- 
ponent lu obtained in both RG schems are given in the 
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second column of Tab. |. For the uncertainty of the four 
loop results we simply take the difference between the 
four loop and the three loop result. This is suggested by 
the behavior of u) in succesive numbers of loops shown 
in Fig. |l} Although both RG schemes lead to compa- 
rable values for to, the convergence of the values in the 
massive scheme is much faster (however the accidentally 
very small error for to derived from this procedure does 
not hold for other quantities). Note that the result for 
lo combined with the correspondipgibur loop results for 
the asymptotic critical exponentaiSEH confirms the con- 
jectured inequality, — vlo < a < 0, for the random mod- 
cls-jcritical exponents involving the specific heat exponent 



massive RG 



-♦ MC [1 



minimal subtraction RG 



wi = 2e + 3. 70401 1 194 e 3/2 + 11. 30873837 e 2 , 
lo 2 = 0.6729265850 e 1/2 - 1.925509085 e 
- 0.5725251806 e 3/2 - 13.93125952 e 2 . 



(7) 



From naively adding the successive perturbational con- 
tributions one observes that already in three loop ap- 
proximation (~ e) 102 becomes negative and therefore no 
stable fixed point exists in strict -^-expansion. Even 
the resummation procedures we applied above, do not 
change this picturecll. This can be considered as indirect 
evidence that the y^-expansion is not Borel summable, 
as may be expected from Ref£3. A physical reason might 
be the existence of the Griffith singularities caused-by the 
zeros of the partition function of the purepgystemES. The 
fixed d approach, both within the massive£3 and minimal 
subtraction! 3 }!!-!! schemes, seems to be the only reliable 
way to study critical behaviour of the model by means of 
RG technique. 
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number of loops 

FIG. 1. Correction-to-scaling exponent lo of the dilute 3d 
Ising model in increasing number of loops. Open square, with 
error bar shows the region of accuracy of the MC dataEI; full 
squares: our values in the minimal subtraction RG scheme; 
full diamonds: our values in the massive RG scheme. 



As it was noted above five loop results for the min- 
imal subtraction scheme are available!!!! In particu- 
lar applying the resumation schemcE3 to the pure Ising 
model case, v — 0, we get the following values for lo 
in increasing number of loops starting from two loop: 
lo = 0.566; 0.852; 0.756; 0.791. This leads to an improve- 
ment in accuracy of the previously calculated d = 3 five 
loop valueE-3 since lo — 0.791 ± 0.036, the value and its 
uncertainty now is comparabler-tp the six loop calcula- 
tion within the massive schemes (see the third column 
of Tab. @). 

The degeneracy of the dilute Ising model /3pfaiactions on 
the one loop level leads to the y^i-expansionU'Ea. For tiie 
critical exponents this expansion is knowp_up to 0(e 2 )cB. 
Starting from the five loop results of RefJij in the replica 
limit we get the following expansionsE3 for the eigenval- 
ues lo\ and lo 2 of the stability matrix (^|) in the fixed point 
u* ^ 0, v* ^ 0: 
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